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Abstract. 

Recent proposed "loophole-free" Bell tests are discussed in the light of classical 
models for the relevant features of optical parametric amplification and homodyne 
detection. The Bell tests themselves are uncontroversial: there are no obvious 
loopholes that might cause bias and hence, if the world does, after all, obey local 
realism, no violation of a Bell inequality will be observed. Interest centres around 
the question of whether or not the proposed criterion for "non-classical" light is 
vàlid. If it is not, then the experiments will fail in their initial concept, since both 
quantum theorists and local realists will agree that we are seeing a purely classical 
effect. The Bell test, though, is not the only criterion by which the quantum- 
mechanical and local realist models can be judged. It is suggested that the 
quantum-mechanical models given in the proposals will also fail in their detailed 
prcdictions. If the experiments are extended by including a range of paramctcr 
vàlues and by analysing, in addition to the proposed digitised voltage differences, 
the raw voltages, the models can be compared in their overall performance and 
plausibility. 
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1. Introduction 

No test of Bell's inequalities EJ to date has been free of "loopholes" . This means 
that, despite the high levels of statistical significance frequently achieved, violations 
of the inequalities could be the effects of experimental bias of one kind or another, not 
evidence for the presence of quantum entanglement. Recent proposed experiments by 
García-Patrón et al |3] and Nha and Carmichael -4 show promise of being genuincly 
free from such problems. If the world in fact obeys local realism, they should noi, 
therefore, infringe any Bell inequality. 




Figure 1. Proposed experimental set- up. In the current text, phase shifts and 
cj> are renamed 8a and 9b- For further explanation, see main text. (Reprinted 
with permission from García-Patrón et al., Phys. Rev. Lett. 93, 130409 (2004). 
Copyright (2004) by the American Physical Society. ) 

The current article discusses a classical model that should be able, once all 
relevant details are known, to explain the results, accounting not only for the failurc 
to infringe the selected Bell test but also for other failures in the detailed predictions. 
It depends on the classical theory for homodyne detection (re-derived here from 
first principies) and the known behaviour of (degenerate-mode) optical parametric 
amplifiers (OPA) 0. 

As far as the "loophole-free" status of the proposed experiments is concerned, 
there would appear to be no problem. A dimculty that seems likely to arise, though, 
is that theorists may not agree that the test beams used were in fact "non-classical" , 
so the failure to infringe a Bell inequality will not in itself be interpreted as showing a 
failure of quantum mechanics§. The criterion to be used to establish the non-classical 
nature of the light is the observation of negative vàlues of the Wigner density, and 
there is reason to think that, even if the Standard method of estimation seems to show 
that these are achieved, the method may be in error. Wigner density is, in any event, 
irrelevant to our model. Far from being, as suggested by García-Patrón and others, 
the "hidden variable" needed, it plays no part whatsoever. 

§ The predicted violation is in any case small, so failure may be put down to other "experimental 
imperfections" . 
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Regardless of the outcome of the Bell tests, and whether or not the light is declared 
to be non-classical, there are features of the experiments that can uscfully be cxploited 
to compare the strengths of quantum mechanics versus (updated) classical theory as 
viable models. The two theories approach the situation from very different angles. 
Classical theory traces the causal relationships between phenomena, starting with the 
simplest assumption and building in random factors later where necessary. Quantum 
mechanics starts with models of complete ensembles, all random factors included. 
This, it is argued, is inappropriate, since two features of the proposed experiments 
demand that we consider the behaviour of individual events, not whole ensembles: 
the process of homodyne detection itself, and the Bell test. 

2. The proposed experiments 

The experimental set- up proposed by García-Patrón et al. is shown in Fig. 1, that of 
Nha and Carmichael being similar. In the words of the García-Patrón et al. proposal: 

The source (controlled by Sophie) is based on a master làser beam, which 
is converted into second harmònic in a nonlinear crystal (SHG). After 
spectral filtering (F), the second harmònic beam pumps an optical parametric 
amplificr (OPA) which generates two-mode squeezed vacuum in modes A and 
B. Singlc photons are conditionally subtracted from modes A and B with the 
use of the beam splittcrs BSa and BSb and single-photon detectors PD^ and 
PDb- Alice (Bob) measures a quadrature of mode A (B) using a balanced 
homodyne detector that consists of a balanced beam splitter BS3 (BS4) and 
a pair of highly-cfficient photodiodes. The local oscillators LOa and LOb 
are extracted from the làser beam by means of two additional beam splitters 
BSi and BS 2 . 

The classical description, working from the same figure, is just a little different. 
Quantum theoretical terms such as "squeezed vacuum" and "quadrature" |j are not 
used since they are not appropriate to the model and would cause confusion. The 
description might run as follows: 

The master làser beam (which is, incidentally, pulsed) is frequency-doubled in 
the crystal SHG. After filtering to remove the original frequeney, the beam is used to 
pump the OPA, a resonant cavity containing a nonlinear crystal cut so as to produce 
degenerate parametric down-conversion of the input. The output comprises pairs of 
classical wave pulses at half the input frequeney, i.e. at the original làser frequeney. 
The selection of pairs for analysis is done by splitting each output at an unbalanced 
beamsplitter (BSa or BSb), the smaller parts going to sensitive detectors PD^ or 
PDb- Only if there are detections at both PD^ and PDb is the corresponding 
homodyne detection included in the analysis. The larger parts proceed to balanced 
homodyne detectors, i.e. ones in which the intensities of local oscillator and test inputs 
are approximately equal. The source of the local oscillators LOa and LO b is the same 
làser that stimulated, after frequeney doubling, the production of the test beams. 

|| The usual model for the elèctric field as the sum of two orthogonal quadratures is appropriate 
where there is no base-line for the phase but not, as here, where all phases concerned are defined 
and measured relativc to a definite base, namely that of the master làser. As will be seen, it it phase 
differences of 180° , not 90° , that feature in the classical model. 
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3. Homodyne detection 

In (balanced) homodyne detection, the test beam is mixed at a beamsplitter with a 
local oscillator beam of the same frequency and the two outputs sent to photodetectors 
that produce voltage readings for every input pulse. In the proposed "loophole-free" 
Bell tests the difïerence between the two voltages will be converted into a digital signal 
by counting all positive vàlues as +1, all negative as -1. 



E 




E t 

Figure 2. Inputs and outputs at the beamsplitter in a homodyne detector. 

is the local oscillator beam, E the test beam, Et and E r the transmitted and 

reflected beams respectively. 

Assuming the inputs are all classical waves of the same frequency and there are 
no losses, it can be shown (see below) that the difference between the intensities of the 
two output beams is proportional to the product of the two input intensities multiplied 
by sin8, where 9 is the phase difference between the test beam and local oscillator. If 
voltages are proportional to intensities then it follows that the voltage difference will 
be proportional to sin#. Whcn digitised, this transforms to a step function, taking 
the value —1 for — n < 9 < and +1 for < 9 < n. (The function is not well defincd 
for integral múltiples of tt.) 

3.1. Classical derivation of the predicted voltage difference 

Assume the test and local oscillator signals have the same frequency, u>, the time- 
dependent part of the test signal being modelled by é 1 ^ , where (ignoring a constant 
phase offset^f) <fi — uit is the phase angle, and the local oscillator phase and test beam 
phascs diffcr by 9. [Note that although complex notation is used here, only the real 
part has meaning: this is an ordinary wave equation, not a quantum-mechanical "wave 
function" . To allay any doubts on this score, the derivation is partially repeated with 
no complex notation in the Appendix.] 

Let the elèctric fields of the test signal, local oscillator and reflected and 
transmitted signals from the beamsplitter have amplitudes E, El, E r and E t 
respectively, as shown in Fig. [21 Then, after allowance for phase delays of 7r/2 at 
each refiection and assuming no losses, we have 

E r = J_(£ e '(0+V2) + ELe *(4>+e)) (!) 

% The phase offset depends on the difference in optical path lcngths, which will not in practice be 
cxactly constant due to thermal oscillations. If a complete model is ever constructed, this should 
therefore be a parameter. 
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and 
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The intensity of the reflected beam is therefore 
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Similarly, it can be shown that the intensity of the transmitted beam is 

1 
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If the voltages registered by the photodetectors are proportional to the intensities, it 
follows that the difference in voltage is proportional to 2£ , £xsin#. When digitised, 
this translates to the step function mentioned above. The probabilities for the two 
possible outcomes are, as shown in Fig. |3 



and 



P+ 



1 for -7T < 6 < 
for < 9 < vr 



for -vr < 9 < 

1 for < 6 < 7T 



(5) 



(6) 



Note that the probabilities are undefincd for integral múltiples of ir. In practice it 
would be reasonable to assume that, due to the presence of noise, all the vàlues were 
0.5, but for the present purposes the integral vàlues will simply be ignored. 
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Figure 3. Probabilities of '+' and '-' outcomes versus phase difference, 
using digitised difference voltages from a perfect, noise-free, balanced homodyne 
detector. 
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4. Application to the proposed Bell tests 

If the freqüències and phases of both test beams and both local oscillators were all 
identical apart from the applied phase shifts, the experiments would be expected to 
produce step function relationships between counts and applied shifts both for the 
individual (singles) counts and for the coincidences. 

It may safely be assumed that this is not what is observed. It would have shown up 
in the preliminary trials on the singles counts (see ref. [7]), which would have followed 
something suggestive of the bàsic predicted step function as the local oscillator phase 
shift was varied. What is observed in practice is more likely to be similar to the 
results obtained by Breitenbach et al. [5] . Their Fig. 6a, reproduced here as Fig. 
shows a distribution of photocurrents that is clustered around zero, for 9 taking integer 
múltiples of 7r, but is scattered fairly equally among positive and negative vàlues in 
between. 




Time [rus] 

Figure 4. A typical scatter of "noise amplitude" (related to voltagc difference) 
as phase 8 is scanned over time. Minimum scatter occurs for 6 taking integral 
múltiples of 7T. (Reprinted from G. Breitenbach et al , J. Opt. Soc. Ara. B 12 
2304 (1995).) 

When digitised, the distribution would reduce to two straight horizontal lincs, 
showing that for each choice of 9 there is an equal chance of a '+' or a '— ' result. As 
in any other Bell test setup, though, the absence of variations in the singles counts does 
not necessarily mean there is no variation in the coincidence rates. As explained in 
the next section, however, the coincidence curves are not the zig-zag ones of Standard 
classical theory. These would be expected if we had full "rotational invariance" + . If 
the ideas presented here are correct, we have instead, in the language of an article 
by the present author jü], only binary rotational invariance. Breitenbach's scatter 
of photocurrent differences is seen as evidence that the relative phase can (under 
perfect conditions) take just one of two vàlues, or 7r. The scatter is formcd from a 
superposition of two sets of points, corresponding to two sine curves that are out of 
phase, together with a considerable amount of noise. 

+ "Rotational invariance" means the hidden variable takes all possible vàlues with equal probability. 
In the current context, if the experiment does indeed produce high visibility coincidence curves, 
the hidden variable responsible will be the common phase difference between test beams and local 
oscillators before application of the phase shifts ( "detector settings" ) 9a and 8b ■ It is argued that in 
the proposed experiment there will be at best approximate rotational invariance, if there is appreciable 
variation in the (again common) frequency. 
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This interpretation accords well with morè comprehensive results of thc 
experiment as reported elsewhere ^ü] . When part of the initial làser beam is used as 
"seed" to the OPA, judicious adjustments of the phase can produce "bright squeezed 
light" and a scatter with alternately positive and negative points. The presence of thc 
seed causes selection of one particular phase set. 

The two "phase sets" arise from the way in which the pulses are produced, 
which involves, after the frequency doubling, the degenerate case of parametric down- 
conversion, the latter producing pulses that are (in contrast to the general case of 
conjugate freqüències) of exactly equal frequency. Consider an initial pump làser of 
frequency lo. In the proposed experiment, this will be doubled in the crystal SHG to 
2ui then halved in the OPA back to w. At the frequency doubling stage, one làser input 
wave peak gives rise to two output ones. Assuming that there are causal mechanisms 
involved, it seems inevitable that every other wave peak of the output will be exactly 
in phase with the input. When we now use this input to produce a down-conversion, 
the outputs will be in phase either with the even or with the odd peaks of the input, 
which will make them either in phase or exactly out of phase with the original làser. 
[The matter can alternatively be approached mathematically, as per ref. where it 
is trcatcd as resonance in a nonlinear situation in which there are two solutions.] 

We thus have two classes of output, differing in phase by ir. If we define the 
random variable a to be for one class, tt for the other, this will clearly be an 
important "hidden variable" of the system. 

The existence of two classes of output of exactly equal frequency and exactly 
opposite phase may well be a feature common to a number of different experiments 
employing degenerate parametric down-conversion sources. One example is discussed 
in ref |ü], namely the Bell test experiment conducted by Weihs et al ^J, but there 
are many more and further examples, not all concerned with Bell tests, will be given 
in later papers. Accepted theory is handicapped by various pre-conceptions. In some 
cases, for example ref. [E], the assumption is made that even in the degenerate case 
the output pair have conjugate, not identical, freqüències. (If this is the case in the 
proposed experiment, though, it will severely reduce the visibility of any coincidence 
curve observed when the experimental beam is mixed back with the source làser in 
the homodyne detector.) In other cases the use of the Standard model in terms of 
orthogonal quadratures leads to neglect of more appropriate models. 

As regards the possibility of any difference in frequency between the test beam 
and the master làser, the preliminary experiments for the García-Patrón proposal 7 , 
using just one output beam may already be sufHcient to show that the interference 
is stronger than would then be the case. It is known that the source làser has quite 
a broad band width, i.e. that lüq is not constant. Though it is likely that it is only 
part of the pump spectrum that induces a down-conversion, so that the band width 
of the test beam may be considerably narrower than that of the pump, it too is non 
zero. It follows that agreement of frequency between this and the test beam must be 
because we are always dealing, in the degenerate case, with exact frequency doubling 
and halving. 

5. A classical model of the proposed Bell test 

In the proposed Bell test of García-Patrón et al. , positive voltage differences will be 
treated as +1, negative as -1. Applying this version of homodyne detection to both 
branches of the experiment, the CHSH test (—2 < S < 2) will then be applied to 
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the coincidence counts. Under quantum theory it is expected that, so long as "non- 
classical" beams are employed, the test will be violated. However, since there are no 
obvious loopholes in the actual Bell test (see later), there is no apparent reason in 
our model why local realism should not win: the test should not be violated. In the 
classical view, this prediction is unrelated to any supposed non-classical nature of the 
light. 



5.1. The bàsic local realist model 

If we take the simplest possible case, in which to all intents and purposes all the 
freqüències involved are the same, the hidden variable in the local realist model is 
clearly going to be the phase difference (pt — or tt) between the test signal and the 
local oscillator. If high visibility coincidence curves are seen, it must be because the 
vàlues of a are identical for the A and B beams. Assuming no noise, the bàsic model 
is easily written down. 

From equation ©, the probability of a — 1 outcome on side A is 

, . í 1 for -7r < 6 A - a < , , 

P-(^,«) = | iol0< e A _ a<7T} (7) 

where 9a is the phase shift applied to the local oscillator A, a is the hidden variable 
and all angles are reduced modulo 2ir. Similarly, the probability of a +1 outcome is 

, . í for — 7r < 9a - ol < , . 

P+(^.«) = ( l iorO<6 A -a<n, ® 

Assuming equal probability \ for each of the two possible vàlues of a, the Standard 
"local realist" assumption that independent probabilities can be multiplied to give 
coincidence ones leads to a predicted coincidence rate of 

P ++ (6 A ,9 B ) = ±p + (9 A ,0)p+(9 B ,0) 

+ -p + (9a,tt)p+(9b,k), (9) 

with similar expressions for , P |_ and P 

The result for 9a — k/2, for example, is 

for -7T < 6 B < 



r ++ (-/W*) = { 1/2 ior0 <9 B <,. ^ 

For 9a = -tt/2 it is 

P ++ (-7T/2,fe)-| Q for0<fe<7r> (11) 

Note that, as illustrated in Fig. the coincidence probabilities cannot, in this 
bàsic model, be expressed as functions of the difference in detector settings, 9 B — 9a- 
This failure, marking a significant deviation from the quantum mechanical prediction, 
is an inevitable consequence of the fact that we have (as mentioned earlier) only binary, 
not full, rotational invariance. 
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Figure 5. Predicted coincidence curves for the ideal experiment, (a) and (b) 
illustrate the settings most likely to be chosen in practice, giving the strongest 
correlations. 9 a is fixed at 7r/2 or —n/2 while 8g varies. In theory, any value 
of 8a between and tt would give the same curve as (a), any between — ir and 
the same as (b). An example is shown in (c), where 8a is 7r/4 but the curve is 
identical to (a). We do not have rotational invariance: the curve is not a function 
of 8 B -8 A . 

5.2. Fine-tuning the model 

Many practical considerations mean that the final local realist prediction will probably 
not look mudi like the above step function. It may not even be quite periòdic. The 
main logical diffcrence is that, despite all that has been said so far, the actual variable 
that is set for the local oscillators is not directly the phase shift but the path length, 
and, since the frequency is likely to vary slightly from one signal to the next (though 
always keeping the same as that of the pump làser), the actual phase difference between 
test and local oscillator beams will depend on the path length difference and on the 
frequency. In a complete model, therefore, the important parameters will be path 
length and frequency, with phase derived from these. 

If frequency variations are sufhciently large, the situation may approach one of 
rotational invariance (RI), but it seems on the face of it unlikely that this can be 
achieved without loss of correlation. If we do have RI, perhaps produced artificially 
by introducing random phase variations into the OPA pump beam, the model becomes 
the Standard realist one in which the predicted quantum correlation varies linearly with 
difference in phase settings, but it is more likely that what will be found is curves that 
are not independent of the choice of individual phase setting. They will be basically 
the predicted step functions but converted to curves as the result of the addition of 
noise. 

It is essential to know the actual experimental conditions. Several relevant factors 
can be discovered by careful analysis of the variations in the raw voltages in each 
homodyne detection system. If noise is low, the presence of the two phase sets, and 
whether or not they are equally represented, should become apparent. 

All this complcxity, though, has no bearing on the bàsic fact of the existence of 
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p++(7c/2,e B ) 




Figure 6. Likely appearance of coincidence curves in a real experiment with 
moderate noise. 



a hidden variable model and the conseqüent prediction that the CHSH Bell test will 
not be violated. 



5.3. The role of the "event-ready" detectors 

In the quantum-mechanical theory, the expectation of violation of the Bell test all 
hinges on the production of "non-classical" light. The light directly output from the 
crystal OPA is assumed to be Gaussian, i.e. it takes the form of pulses of light that 
have a Gaussian intensity profile and also, as a result of Fourier theory, a Gaussian 
spectrum. When this is passed through an unbalanced beamsplitter (BSa or BSs) 
and a "single photon" detected by an "event-ready" detector, the theory says that 
the subtraction of one photon leaves the remaining beam "non-Gaussian" . Although 
there is mention here of single photons, the theory is concerned with the ensemble 
properties of the complete beams, not with the individual properties of its constituent 
pulses. 

In the local realist (classical) model considered here, the shapes of the spectra are 
not relevant except insofar as a narrow band width is desirable for the demonstration 
of dramàtic correlations. The "event-ready detectors" (PD^ and PD# in Fig.QJ play, 
instcad, the important role of selecting for analysis only the strongest down-converted 
output pairs, it being assumed that the properties of the transmitted and reflected 
light at the unbalanced beamsplitters are identical apart from their amplitudes. It is 
likely that those detected signals that are coincident with each other will be genuine 
"degenerate" ones, i.e. of exactly equal frequency, quasi-resonant with the pump làser. 
The unbalanced beamsplitters and the detectors PD^ and PD b need to be set so that 
the intensity of the detected part is sumcicnt to be above the minimum for detection 
but low enough to ensure that all but the strongest pulses are ignored. 

In neither theory are the event-ready detectors really needed in their "Bell test" 
role of cnsuring a fair test (sec below), since the homodyne detectors are almost 100% 
cfBcicnt. 



6. Validity of the proposed Bell test 

Coincidences between the digitised voltage differences will be used in the CHSH 
Bell test ^3 EL but avoiding the "post-selection" that has, since Aspect's second 
experiment |15j . become customary. The García-Patrón et al. proposal is to use 
event-ready detectors, as recommended by Bell himself for use in real experiments 
|16| . None of the usual loopholes ^Ij are expected to be applicable: 

(i) With the use of the event-ready detectors, non-detections are of little concern. 
The detectors (PD^ and PDb in Fig. act to define the sample to be 
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analysed, and the fact that they do so quite independently of whether or 
not any member of the sample is then also detected in coincidence at the 
homodyne detectors ensures that no bias is introduced here. The estimate of 
"quantum correlation" * to be used in calculating the CHSH test statistic is 
E = (N ++ ab + N--ab — N + -ab — N- + ab)/Nab, where the N's are coincidence 
counts and the subscripts are self-explanatory. This contrasts with the usual 
method, in which the denominator used is not Nab but the sum of observed 

coincidences, N ++ ab + N ab + ab + N ^ab- The use of the latter can 

readily be shown to introduce bias unless it can be assumed that the sample of 
detected pairs is a fair one. That such an assumption fails in some plausible local 
realist models has been known since 1970 or earlier |14II18| . 

(ii) There is no possibility of synchronisation problems Jü] , since a pulsed source is 
used. 

(iii) No "accidentals" will be subtracted |17|. 

(iv) The "locality" loophole can be closed by using long paths and a random system 
for choosing the "detector settings" (local oscillator phase shifts) during the 
propagation of the signals. 

The system is almost certainly not going to be "rotationally invariant" (not all 
phase differences will be equally likely), but this will not invalidate the Bell test. It 
may, however, be important in another way. It is likely that high visibilities will be 
observed in the coincidence curves (i.e. high vàlues of (max - min)/(max + min) in 
plots of coincidence rate against difference in phase shift), leading to the impression 
that the Bell test ought to be violated. These visibilities, though, will depend on the 
absolute vàlues of the individual settings. High ones will be balanced by low, with the 
net effect that violation does not in fact happen. 

7. Validity and significance of negative estimates for Wigner densities 

In the quantum mechanical theory discussed in the loophole-free Bell test proposals 
and in other recent papers [HI d, P ar t of the evidence that is put forward as indicating 
that negative Wigner densities are likely to be obtained consists in the observation 
that, when 9 is varied randomly, the distribution of observed voltage differences shows 
a double peak (see Fig. [7J. There is a tendency to observe roughly equal numbers of 
+ and - results but relatively few near zero. The fact that the relationship depcnds 
on the sine of 8 is, however, sufhcient to explain why this should be so. 

To illustrate, let us consider the following. The sine of an angle is between and 
0.5 whenever the angle is between and ir/6. It is between 0.5 and 1.0 when the angle 
is between ir/6 and ir/2. Since the second range of angles is twice the first yet the 
range of the vàlues of the sine is the same, it follows that if all angles in the range to 
7r/2 are selected equally often there will be twice as many sine vàlues seen above 0.5 
as below. The same holds when random angles between n/2 and ir are chosen, whilst 
for vàlues between —tt and we find a symmetrical result for negative vàlues. 

When allowance is made for the addition of noise, the production of a distribution 
such as that of Fig. [3 for the average when angles are sampled uniformly comes as no 

* In the derivation of Bell inequalities such as the CHSH inequality, the statistic required is simply 
the expectation value of the product of the outcomes. This is commonly referred to in this context as 
the "quantum correlation" . Only when there are no null outcomes, possible outcomes being restricted 
to just +1 or —1, does it coincide with ordinary statistical correlation. 



Homodyne detection and optical parametric amplification 



12 



x 



-2 







2 








0.2 



0.4 



pr(X) 



Figure 7. Observed distribution of voltage differences, X, using homodyne 
detection in an experiment similar to the proposed Bell test and averaging over a 
large number of applied phase shifts. (Based on Fig. J^a of A. I. Lvovsky et al., 
Phys. Rev. Lett. 87, 050402 (2001).) 



surprise. Clearly, as the experimenters themselves recognise, the dip is not in itself 
sufficient to prové the non-classical state of the light. For this, direct measurement 
of the Wigner density is required, but there is a problem here. No actual direct 
measurement is possible, so it has to be estimated, and the method proposed is 
the Radon transformation [20]. It is claimed that in other experiments Wigner 
densities calculated either by this procedure or by a "more precise maximum-liklihood 
reconstruction technique" [22] have shown negative regions, but perhaps the methods 
should be checked for validity? 

In any event, as already explained, the natural hidden variable relevant to the 
proposed experiment is the phase of the individual pulse, not any statistical property 
of the whole ensemble. Indeed, the use of Wigner density as a substitute for hidden 
variables was never originally intended and has no theoretical basis. In Bell's much- 
quoted paper on the subject (Ch. 21 of ref. |23)> the "hidden variables" remain, 
as ever, parameters such as position and momentum that are specific to individual 
partides. The role of the negative Wigner density is merely to provide, in certain 
rather special circumstances, and alternative test for nonlocality in that, if negative 
vàlues are found, then real local hidden variables cannot exist. 

8. Suggestions for extending the experiment 

The bàsic set-up would seem to present an ideal opportunity for investigation of some 
key aspects of the quantum and classical models, as well as the operation of the Bell 
test "detection loopholc" . 

(i) The operation of the detection loophole could be illustrated if, instead of 
using the digitised difference voltages of the homodyne detectors, the two separate 
voltages are passed through discriminators. The latter operate by applying a 
threshold voltage that can be set by the experimenter and counting those pulses 
that exceed it. These can be used in a conventional CHSH Bell test, i.e. using total 
observed coincidence count as denominator in the estimated quantum correlations 
E. The model that has been known since Pearle's time (1970) prediets that, as the 
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threshold voltage used in the discriminators is increased and hence the number 
of registered events decreased (interpreted in quantum theory as the detection 
efficiency being decreased), the CHSH test statistic S, if calculated using estimates 
E = (N ++ + 7V__ - jV+_ - N- + )/(N ++ + 7V__ + iV+_ + N-+), will increase. If 
noise levels are low, it may well exceed the Bell limit of 2. 

Such an experiment has, in a sense, already been performed by Babichev et al |22| 
and yielded the expected results. If it is accepted that their source (two outputs 
from a beamsplitter) would have been in an entangled state, then their Fig. 4b 
clearly demonstrates that high detector thresholds (i.e. low detector efficiencies) 
can lead to violations of the Standard form of the CHSH test. 

(ii) The existence of the two phase sets is, in point of fact, well known when 
OPA's are operated above threshold The resonance into one or other of the 
sets then becomes stable. The fact that that two sets are also responsible for many 
of the observations below threshold could be further investigated if either the raw 
voltages or the undigitised difference voltages are analysed. So long as the noise 
level is low, the existence of the two superposed curves, one for a = and the 
other for a — tt, should be apparent. It would be interesting to investigate how 
the pattern changed as optical path lengths were varied. Breitenbach's pattern 
might be hard to reproduce using long path lengths, where exact equality is 
needed unless the light is monochromatic. 

(iii) Comparison of overall performance: If the primary goal of the experimenter 
is clearly set out to be the comparison of the performance of the two rival models, 
rather than merely the conduct of a Bell test, further ideas for modifying the set- 
up will doubtless emerge when the first experiments have been done. Many of 
the predictions of the quantum-mechanical model have already appeared in print 
PJE]- The first stage in comparing models should probably be, therefore, to 
conduct supplcmentary experiments so as to establish the relevant parameters 
of the full classical model and hence make equivalent empirical predictions. 
It is possible, though, that qualitative predictions alone will be sufficicnt to 
demonstrate superiority one way or the other. 

9. Conclusion 

The proposed experiments would, if the "non-classicality" of the light could be 
demonstrated satisfactorily, provide a definite answer one way or the other regarding 
the reality of quantum entanglemcnt. Thcy could usefully be extended to includc 
empirical investigations into the operation of the Bell test detection loopholc. Pcrhaps 
more importantly, though, they present valuable opportunities to compare the 
performance of the two theories in both their total predictive power and their 
comprehensibility. Are parameters such as "Wigner density" and "degree of squeezing" 
really the relevant ones, or would we gain more insight into the situation by talking only 
of freqüències, phases and intensities? Parameters such as the detection efficiency and 
the transmittance of the beamsplitters will undoubtedly affect the results, but do they 
do this in the way the quantum theoretical model suggests? It will take considerably 
more than just the minimum runs needed for the Bell test if we are to find the answers. 

The detailed predictions of the classical model cannot be given until the full 
facts of the experimental set-up and the performance of the various parts are known, 
but it gives, in any event, a simple explanation of the double-peaked nature of the 
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distribution of voltage differences. The peaks arise naturally from the way in which 
homodyne detection works, and the quantum theoretical idea that they are one of the 
indications of a non-classical beam or of negative Wigner density would not appear to 
be justifiablc. The idea that a classical beam can become non-classical by the act of 
"subtracting a photon" is, equally, of doubtful validity. The experimental role of the 
subtraction and detection of part of each beam is to aid the selection for coincidence 
analysis of those pulses that are likely to be most strongly correlated. 
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Appendix. Alternative classical derivation of the homodyne detection 
formula 

A derivation is given here that does not involve complex numbers and hence confirms 
that the equations in the text are ordinary wave equations, not quantum- mechanical 
"wave functions" . 

The relationship between intensity difference and the local oscillator phase can 
be checked as follows: 

Assume the two input beams are 

Experimental beam: E cos (f> , where, as before, cf> = ujt 
Local oscillator: El cos(<^> + 9) 

Then the output beams, assuming a 50-50 beamsplitter and no losses, can be written 

Rcflccted beam: 

E r = -j=(E cos(> + n/2) + E L cos(> + 9)) 

= -=(-Esm.^ + E L cos(^> + 9)) (A.l) 

Transmitted beam: 

E t = —j= (Ecosct> + E L cos((f> + + Tr/2)) 

= — (£cosç6- E L sm((f> + 6)). (A.2) 

Let us define a (constant) angle tp such that tan í/j = E/El, making Ecostp = 
El siiiip. 

Consider the case when 9 = 0. We have 

E r = —=(E/ sin·0)(— smíf) sin0 + cos-0 cos (f>) 
v2 

= ^=(E/smi/))cos(ilj + <f)), (A.3) 
v2 
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so that amplitude is proportional to -^=E/ simp and intensity to \E 2 j sin 2 ip. 

The intensity of the transmitted beam can be found similarly, and turns out to 
be also proportional to \E 2 / sin 2 tjj, so that the voltage difference from the homodyne 
detector is therefore to be zero. Likewise, a zero difference is found for 8 — w, but 
other vàlues of 9 produce more interesting results. 

For example, for 9 — n/2 we find 

E r = —=(Ej sin ?/>)(— s'mi/j sin (f> — cosip sm(f>) 
V2 

= —=(Ej s'mijj) sm<j)(— smip — cosip), (A.4) 
v2 

Et = —p (E I sin tp) (sin 0> cos <f> — cos i\) cos <f>) 
V2 

= —=(E/smip) coscj)(smip — cos-0). (A.5) 
The difference in intensities is therefore proportional to 

Difference = -(i? 2 / sin 2 ?/')[(— sin i/j — cos ip) 2 
— (sinV' — cosí/>) ] 
= - (E 2 / sin 2 sin ^ cos ip) 

= —E 2 cos-0/ sin?/; 
= 2E 2 E L /E 

= 2EE L . (A.6) 

This is consistent with the result obtained by the method in the main text, so it 
seems safe to accept that as being correct. 
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